A Dark-Matter Spike at the Galactic Center? 
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The past growth of the central black hole (BH) might have enhanced the density of cold dark matter 
halo particles at the Galactic center. We compute this effect in realistic growth models of the present 
(2 - 3) x lO 6 Af BH from a low -mass seed BH, with special attention to dynamical modeling in a 
realistic galaxy environment with merger and orbital decay of a seed BH formed generally outside 
the exact center of the halo. An intriguing "very-dense spike" of dark matter has been claimed in 
models of Gondolo and Silk with density high enough to contradict with experimental upper bounds 
of neutralino annihilation radiation. This "spike" disappears completely or is greatly weakened when 
we include important dynamical processes neglected in their idealized/restrictive picture with cold 
particles surrounding an at-the-center zero-seed adiabaticly-growing BH. For the seed BH to spiral 
in and settle to the center within a Hubble time by dynamical friction, the seed mass must be at 
least a significant fraction of the present BH. Any subsequent at-the-center growth of the BH and 
steepening of the central Keplerian potential well can squeeze the halo density distribution only 
mildly, whether the squeezing happens adiabatically or instantaneously. 



PACS numbers: 95.35.+d; 14.80.Ly; 95.55.Vj 
I. INTRODUCTION 

Weakly interacting massive particles (WIMPs) are 
among the leading candidates for the dark matter in 
galactic halos. Such particles arise in extensions to the 
standard model of particle physics. The most widely 
studied example is the neutralino, plausibly the lightest 
superpartner in supersymmetric versions of the standard 
model. The cosmological abundance of WIMPs is natu- 
rally of order the critical density of the Universe, and so 
WIMPs make natural dark-matter candidates. 

The search for WIMPs in the Milky Way halo has been 
a major endeavor both theoretically and experimentally 
in the last two decades (for comprehensive reviews see 
Refs. UK). The main efforts so far have been to de- 
tect WIMPs directly in low-background laboratory ex- 
periments and indirectly via observation of the energetic 
neutrinos produced by annihilation of WIMPs that have 
accreted in the Sun and/or Earth. Another possibility 
is to search for exotic cosmic rays, such as high-energy 
gamma rays, positrons, antiprotons, or neutrinos pro- 
duced by WIMP pair annihilations in the Galactic halo. 
In particular, the flux of gamma rays in a given direction 
is proportional to the line-of-sight integral of the square 
of the WIMP density. Since the dark-matter density is 
expected to be largest towards the Galactic center, the 
flux of such exotic gamma rays should be highest in that 
direction. 

We now know that black holes of mass of 10 6 — 10 8 M 
reside in most galactic centers (see, e.g., ||). In partic- 
ular, the arguments for a 2.6 x 10 6 M Q black hole at 
the center of the Milky Way are at this point almost in- 
controvertible (see, e.g., Refs. [Q) given the rise of the 
velocity (radial and proper motion) dispersion within 1 
pc of the Galactic center. It has long been argued that 



if a black hole grows inside a given population of stars, 
the density of stars should be enhanced in the vicinity of 
the black hole after its appearance 0-B. If the WIMP 
density is likewise enhanced at the Galactic center, then 
the source of exotic cosmic rays from WIMP annihilation 
therein should be increased. Several authors investigated 
the enhancement in the gamma-ray flux for the dark- 
matter-density spike, p(r) cx r^ 15 as the radius r — » 0, 
that would be induced by the adiabatic growth of a black 
hole in an initially isothermal dark-matter profile Q. 

More recently, Gondolo and Silk [Bl (hereafter GS) in- 
vestigated the dark-matter spike that arises if a black hole 
grows adiabatically at the center of a dark-matter halo 
that initially has a singular power-law cusp, p(r) oc r 1 
with < 7 < 2, as suggested by numerical simulations 
]lO| , p| . They found that the black-hole growth gives rise 
to a dark-matter spike p(r) oc r~ A with 2.25 < A < 2.5, 
in agreement with earlier scaling results found in Ref . . 
In this case, the cosmic-ray source from WIMP pair anni- 
hilation, proportional to J d 3 r p 2 (r), becomes huge due 
to the diverging contribution from small radii. (The di- 
vergence is cut off by the black-hole horizon or by a max- 
imum density set by WIMP pair annihilation.) The flux 
of neutrinos from neutralino annihilation in this spike ex- 
ceeds current experimental upper bounds in some regions 
of neutralino parameter space. Subsequent work p2 13 
shows that synchrotron emission from the motion of elec- 
trons and positrons produced by neutralino annihilation 
in this spike would also be significant and in most cases 
exceeds current upper bounds. The conclusion of Refs. 
P, ^2|p^| is that either neutralino dark matter or a cuspy 
halo must be ruled out. 

Given these very important implications for dark- 
matter searches, the claim of the possible presence of 
a steep dark-matter spike at the Galactic center and 
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the underlying assumptions should be investigated more 
carefully. Here we argue that although a possibility, the 
enhancement found by GS requires somewhat unusual 
initial conditions for the Galactic halo and for the black 
hole. In particular, we emphasize how crucial it is that 
the black hole grew adiabatically from a tiny initial mass 
to its present-day mass, and that this happened precisely 
at the center of the dark-matter distribution. 

To do so, we consider a few alternatives for the black- 
hole- format ion scenario that we can treat with semi- 
analytic methods and jet captures the main dynamical 
processes in (hierarachical) galaxy formation. We show 
that for the GS spike to arise, the majority of the black- 
hole growth had to occur within the inner 50 pc of the 
dark-matter cusp and on a timescale long compared with 
10 7 years. We also show that the growth of the spike re- 
lies on the existence of a very high density of cold orbits 
within the inner 50 pc. If the black hole is formed away 
from the center and then spirals in, it should gcncrically 



lead to the weaker dark-matter spike, p oc 



--0.5 



with 



possibly even a depletion, rather than an enhancement, 
in the central dark-matter density. 

The outline of the paper is as follows: In Section 0, we 
consider a simple toy model to illustrate the origin of the 
spike found by GS in case the black hole grows adiabat- 
ically at the exact center of the dark-matter system. In 
Section III we perform the computation of the adiabatic 
growth for a class of more realistic dark-matter-halo pro- 
files, stressing a few more subtle points. In Section IV we 
treat instead the opposite case of fast growth for the black 
hole. Section |v| deals with the more complex scenarios 
of off-centered formation for the black hole. While in all 
cases above we include just dark matter and the black 
hole, in Section VI we consider the more realistic sce- 



nario in which other baryonic components (i.e., stars in 
the stellar nucl eus, bulge, and disk) are included as well. 
In Section VII we summarize and comment on current 



theoretical ideas about the origin of the black hole. 



II. WARMUP: HALO DISTRIBUTIONS WITH 
ONLY CIRCULAR ORBITS 

Before proceeding with our numerical calculations, let 
us consider a simple toy model that illustrates the ori- 
gin of the GS spike. In this model, we assume an ini- 
tial self-gravitating dark-matter distribution pi{r) made 
of particles that are all on circular orbits, at the center 
of which a black hole grows adiabatically. The angular 
momentum of each particle remains invariant since the 
central black hole exerts no torque on any dark-matter 
particle. This tells us that TiMiiji) — rfMf(rf), where 
M(r) is the total (i.e., dark-matter plus black-hole) mass 
enclosed within the radius r. Moreover, conservation of 
the dark-matter mass tells us that M- [(rij — Mj?(ry), 
where the superscript h denotes the halo contribution to 
the mass. 



Now suppose the initial dark-matter density profile is 
Pi(ri) oc r~ 7 and let the final density be Pf{rf) oc rJ A . 

Conservation of mass implies that r 3-7 oc r^~ A . In the fi- 
nal state, the mass enclosed within a radius r will be dom- 
inated by the black hole, Mf(rf) ~ Mbh, as r — > 0. Con- 
servation of angular momentum thus gives rf~ 7 oc Tf. 
Putting these two conditions together, we find a final 
profile, p f (r f ) oc rj A , where A = (9 - 2 7 )/(4 - 7) [§,§. 
Note that for < 7 < 2, we get 2.25 < A < 2.5. Hence 
the WIMP pair annihilation rate, which is proportional 
to the density Pf(rf) squared, would have diverged ex- 
cept for the truncation at the horizon of the BH. 

This derivation illustrates several important points, 
some of which will be discussed further for more general 
distributions in subsequent Sections. 



A. Where does the mass in the spike come from? 

The dark-matter particles that make up the spike come 
from the inner part of the dark-matter halo, the part that 
encloses a mass roughly equal to Mbh — 2.6 x 10 6 M Q . 
The addition of a black hole will have negligible effect on 
circular orbits that enclose larger masses. That the spike 
comes from this inner part agrees with intuitions, and is 
suggested by numerical simulation of Ref. § (see their 
Fig. 6) for models with isotropic velocity distributions. 

For the power-law density profile that GS considered 
(and for the NFW profile we consider below), the radius 
that encloses 2.6 x 10 6 M of halo mass is only about 
50-60 pc. Thus, for the spike in the GS picture to form, 
the black hole must have originated within the inner 50- 
60 pc of the Galactic center. This is in some sense sur- 
prising because this length scale is microscopic compared 
with characteristic Galactic distance scales: ~ 3 kpc for 
the bulge and ~ 10 kpc for the disk. Put another way, 
the inner 50 pc encloses only about (1.5 — 2) x 10 8 M Q 
of stellar material, about 5 x 10~ 3 of the total stellar 
mass (probably about (5 - 8) x 10 10 M ) of the galaxy. 
Although the alignment between the black hole and the 
stellar distribution is quite good in our own Galaxy (even 
that is a mystery — the situation is quite different, e.g., 
in M31), the required alignment is between the black 
hole and the center of the dark-matter system. Although 
baryons and dark matter are mixed homogeneously in the 
early Universe, there are numerous processes involved in 
galaxy formation that should segregate the baryons from 
the dark matter. For example, baryons will experience 
radiation drag from the cosmic microwave background 
when they first start to cool; UV radiation from the first 
generation of stars could exert a wind on the baryons; 
in galaxy mergers, baryons shock while dark-matter ha- 
los pass through each other; baryons today appear in 
disks, bulges, and bars, unlike dark matter; they are 
clumped into globular clusters, open clusters, and molec- 
ular clouds. Thus, there are many good reasons to be- 
lieve that the precise center of the baryon distribution 
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should be displaced relative to that of the dark matter, 
at least at some intermediate stages in the formation of 
our Galaxy. Thus, even if the black hole formed at the 
center of a stellar distribution, it would still be some- 
what surprising if, at all times during its formation, that 
coincided to 50 pc with the center of the dark-matter 
distribution. 



B. Cold particles 

The circular-orbit toy model shows us that the parti- 
cles that give rise to the spike are initially very cold; i.e., 
the circular speeds in the initial distribution approach 
zero as r — » 0. (Below, we will show that the same is 
true for more general and realistic velocity distributions.) 
Again, this is somewhat surprising, as cold particles are 
easily subject to heating by a variety of phenomena, es- 
pecially in the environment at the Galactic center. For 
example, the passing of molecular clouds, open, or glob- 
ular clusters would tend to disrupt these cold orbits, as 
would the existence of a bar. Likewise, the mergers and 
sub-halo accretion that certainly played a role in the for- 
mation of the Milky Way would tend to disrupt these cold 
orbits. We cannot be more quantitative without a more 
precise model for galaxy formation. However, these dis- 
ruptive processes must have occurred at some level, and 
they would tend to hinder the growth of the dark-matter 
spike. 

Note that a dark-matter spike that will give rise to a di- 
vergent WIMP annihilation rate can arise even if 7 = 0. 
As explained in Ref. ||, implied above, and discussed 
further below, the existence of a spike is more generally 
due to the existence of a singularity in the density of cold 
orbits at the Galactic center. Such singularities can arise 
even if the density profile is non-singular. If the density 
p at the origin is finite and can be written as an analytic 
function of the Cartesian coordinates x, y, and z, then 
there will be no singularity in the phase-space distribu- 
tion that gives rise to a spike, but otherwise there will 
be a phase-space singularity and a spike after adiabatic 
growth of a black hole H. Thus, consider for example, 
the profiles p(r) oc (a + rj 2 and p(r) cx (a 2 + r 2 ) -1 , both 
of which approach a constant at the origin, and suppose 
that both have isotropic velocity distributions. As ar- 
gued in Ref. || and described below, adiabatic black-hole 
growth in the former gives rise to a spike p cx r~ A with 
A = 2, while the latter gives rise to a spike with A = 1.5. 

C. Initial black- hole mass 

The circular-orbit model also indicates that the spike 
will grow only if the black hole grows almost entirely 
(not just partially) at the black-hole center. Put another 
way, the spike that grows from an initial state with a 
10 6 Mq black hole and a p oc r~ x dark-matter halo to a 



final state with a 2 x 10 6 M Q black hole is not nearly as 
dramatic as that which arises when a black hole is grown 
from nothing to 10 6 Mq. If we start with a black hole of 
significant mass and a dark-matter cusp, then the angular 
momentum of the inner orbits is dominated by the initial 
black-hole mass. Slowly doubling the black-hole mass 
reduces the radius of each of these inner orbits by a factor 
of 2 and thus leads to a factor-of-4 density enhancement. 
In other words, the huge densities of the innermost orbits 
that give rise to a huge annihilation signal come from the 
fact that initially, Mj(r) — ► as r — > and the radius Ti 
of each inner orbit is reduced by a factor ~ Msa/Mi{ri). 
If a black hole exists at the Galactic center from the 
start with a mass Afg H and then grows to a final mass 
Afg H , the initial dark-matter density profile Pi(ri) oc v\ 7 
becomes a density profile with the same exponent and 
normalization enhanced by a factor ~ (M^/M^) 3 ^ 1 . 

D. Timescale for adiabatic growth 

If the initial halo profile is p(r) ~ r _1 , then the en- 
closed mass will be M(r) oc r 2 , the circular velocity will 
be v c (r) oc r 1 / 2 and the period of a circular orbit will be 
T oc r 1 / 2 . The period of an orbit at 50 pc will there- 
fore be ~ 10 7 yr. Thus, the black hole can grow rel- 
atively quickly compared with a cosmological timescale 
and still be considered adiabatically. This provides an 
argument in favor of the GS spike. On the other hand, 
the present-day growth rate of the black hole due to con- 
sumption of stars coming inside its loss cone is at most 
3 x 10 4 M Q /Gyr Jl4| , one order of magnitude smaller 
than what is needed to form the current black hole in 
a Hubble time (mean growth rate of 2 x 10 5 M /Gyr). 
This information combined with the evidence that peak 
AGN activity occurred near redshift z ~ 2 suggests that 
the majority of the black-hole mass was in place at these 
early epochs. Although this leaves plenty of time (few 
Gyr) for the black hole to grow adiabatically, the violent 
processes that are likely required to form the black holes 
that power AGN may well operate on a much more rapid 
timescale. 



E. Weaker spikes at smooth densities 

All of the discussion above should make it clear that if 
the black hole grows adiabatically away from the Galactic 
center, or in a Galactic center with a nonsingular phase- 
space distribution, then the spike should be the weaker 
p oc r~ 15 spike found in Refs. |5|-|^]. Although this will 
produce an enhancement in the flux of annihilation ra- 
diation, it will not produce the huge enhancement that 
arises in a spike p(r) oc r~ A with A > 2. 
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III. ADIABATIC GROWTH OF THE BLACK 
HOLE 

We now consider adiabatic growth of a black hole at the 
center of a spherically symmetric self-gravitating dark- 
matter system with a more general (and realistic) phase- 
space distribution. Because of the two key ingredients 
here — adiabaticity and spherical symmetry — this prob- 
lem can be treated semi-analytically. The slow changes 
in the potential allow us to relate in a simple way the 
dark-matter distribution functions before and after the 
appearance of the black hole. The symmetry, both in 
the initial and final states, implies simple forms for the 
adiabatic invariants. 

The calculation has been discussed in several papers 
in the past, in the context of growth of a black hole in a 
general stellar system |5|-|^], or in the context of growth 
of the central black hole at the Milky Way center || (see 
also Ref. 0). Below we briefly review the calculation, 
both to reproduce the result in GS and to provide some 
additional insight on the origin of the spike. 



A. Semi-analytic treatment 

In a spherically symmetric problem, the adiabatic in- 
variants are the angular momentum and the radial ac- 
tion; i.e^, Lf = Li = L and J r j(Ef, L) — J rt i(Ei, L), 
where [[15| 



J r (E,L) = - 

7T 

_ 1 
7T 



dr v r 



dr \ 2(E - $(r)) 



L 2 



(1) 



and the subscripts i and / refer to the value of the quan- 
tity in the initial and final state. If the growth of the 
black hole is adiabatic, the distribution function f(E,L) 
remains unchanged, so ff(Ef,L) = fi(Ei,L). 

We specify the initial state by its density Pi{r) as a 
function of radius r. We also postulate that it has an 
isotropic velocity-dispersion tensor. If so, then there is 
a one-to-one correspondence between Pi(r) and fi(Ei), 
given by Eddington's formula (see, e.g., |15|), 



f(E) 



1 



V8vr 2 



d 2 



P 



(dp/d^ =Q 



E d^ 2 V$ - E 



, (2) 



and the initial state is fully defined. 

By inverting the expression for the radial action in the 
initial state, we can write Ei as a function of Ef and 
L and hence find ff(E/,L) from fi(Ei,L). To compute 
the radial action in the final state, we actually need first 
a guess for the final potential. We assume that the po- 
tential is dominated by the the black hole in its vicinity, 
while it is roughly unchanged at large galactocentric dis- 
tances, fixing <&/(r) = 3>i(r)-GM B ii/r. Once ff(Ef, Lf) 
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FIG. 1. Enhancement of an NFW dark-matter-density pro- 
file due to the growth of the Galactic black hole at the center 
of the dark-matter system. The initial profile (dotted curve) 
is modified into the solid curve if the growth of the black hole 
is adiabatic, or into the dash-dot curve for a sudden growth. 
Also shown is the case of adiabatic growth for an initial profile 
with all particles on circular orbits (dashed curve). p CO re is the 
maximum WIMP density above which WIMPs are depleted 
by pair annihilations. 

is obtained, the density profile of the system in the fi- 
nal state is recovered by integrating ff over the velocity 
space, or, with a change of variables, over Ef and L: 



Pf( r ) 



4tt 



dE 



,//,, -_/>;/• ,./,,;. ( 3 ) 



Because of the form we choose for <&/(r), we expect Pf(r) 
to be fairly accurate in the central region and far from 
the black hole. A more accurate result could be obtained 
in the transition region by iterating. However, we will 
mainly be interested in the region near the black hole 
where the approximation $/(r) ~ —GMbb/v is fairly 
accurate. It is also the form adopted by GS. 

We apply this procedure numerically in one case. 
Among the general family of spherical double-power-law 
dark-matter profiles (the Zhao models EF 



/ ^0 

Pi(r) = Po [ — 



1 + {R a /a) a 
1 + (r/a) a 



-, (0-y)/a 



(4) 



we focus on the functional form suggested by the N-body 
simulation of Navarro, Frenk, and White pSJ^ , a density 
profile with (a,/3, 7) = (1,3,1) (hereafter the NFW pro- 
file). We fix the values of the Sun's galactocentric dis- 
tance, the local halo density and the length scale in the 
profile, respectively, as Ro = 8 kpc, p — 0.3 GeV cm -3 
and a = 20 kpc (values compatible with available dy- 
namical constraints); this profile is plotted as a dotted 
line in Fig. || 
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We first compute numerically the corresponding distri- 
bution function and then we follow the same steps applied 
in Ref. || for the isothermal sphere. We perform numer- 
ically the inversion between Ei and Ef, computing the 
radial action for the initial and the final state on grids, 
respectively, in the (Ei,L) and (Ef,L) plane (we imple- 
ment adaptive grids in order to achieve everywhere the 
required numerical accuracy for the inversion). Then we 
compute p f (r) performing the double integral in Eq. (||) . 
The result is shown as a solid line in Fig. [|. As can be 
seen, the adiabatic growth of the black hole at the center 
of the dark-matter system induces a sharp increase in the 
dark-matter density close to the black hole (compare the 
solid line with the dotted line which gives the initial pro- 
file). At a galactocentric distance of about 20 pc, there 
is the transition from the 1/r cusp in the NFW profile to 
a much steeper spike p(r) oc r~ A with A = 7/3. This is 
in agreement with the corresponding result (profile scal- 
ing as 1/r 1 with 7 = 1) in GS. In Fig. [j] we also plot, 
as a dashed line, the result for the adiabatic contraction 
if all particles are on circular orbits (the toy model we 
introduced in the previous Section). As can be seen, this 
simplified treatment reproduces quite closely the result 
we obtained in case of an isotropic distribution (within 
a factor of 2, which is less than the uncertainty in the 
overall normalization of the dark-matter profile). 

The density in the spike exceeds 10 s M Q pc~ 3 , which 
is roughly the maximum density if the dark matter is 
made of WIMPs, which we indicate as p C ore in this and 
in the following Figures. The limit comes from the value 
of the pair-annihilation cross section for WIMPs, which 
is fixed, to some extent, by their relic abundance (see 
GS for details). In fact, the cross section for WIMP- 
WIMP elastic scattering should generally be larger than 
the annihilation cross section, so the maximum central 
density derived by GS should actually be a bit lower. 



B. Scaling properties 

The result for the spike slope we have just derived 
numerically could have been anticipated using scaling 
properties, as suggested in Ref. ||. This argument 
can be summarized as follows: Let the potential of the 
initial system increase with radius, for small radii, as 
$i(r) ~ r 2 ~ 7 ; this is the leading term in systems such 
as, e.g., those described by Eq. (Q) with 7 > 0. Then 
the conservation of the radial action implies the scaling 



Ei ~ K 



'(2-7)7(4-7) 



^j. Suppose then that the distribu- 
tion function in the initial state diverges near Ei = $^(0) 
as fi(Ei) — [Ei - <&,(0)] -n . Combining this with the ex- 
pression for the final density profile, Eq. (0), it follows 
that 



Pf( r ) 



a 3 

A = o + n 

2 V 4 -7 



(5) 



the family of profiles introduced above in Eq. (|J). We 
derive n for 1 < a < 2; results for other ranges can be 
obtained analogously. In Eddington's formula, Eq. (0), 
the second term is zero, while the main contribution to 
the first comes from the singularity in the lower extreme 
of integration. After some algebra, one finds that the 
leading terms are given by 



f(E) 



dr 



jcp + Cl (r/a) a + c 2 (r/a) 2 " + . . 
r 3 ^ r 2-7 -[E- $(0)] 



(6) 



where r* = [E— ^(O)] 1 ^ 2 7 ', and the first two coefficients 
in the expansion are (up to a constant), 



c = -27(3 - 7) 

Cl = (a-2)(3- 7 )(/3-7) 

(2a — 7) (3 — 7) + a(a — 7) 
a(3 + a — 7) 



(7) 



To illustrate, consider three cases. For 7 > and inde- 
pendent of a, cq 7^ and 



6 — 7 
2(2-7) 



A = 



9-27 
4-7 



(8) 



recovering the result of Ref. g and GS. In the limit 7^0 
and for I < a < 2, the first term in the expansion van- 
ishes, c = 0, and hence the second becomes the leading 
term. In this case one finds milder singularities, 



.4 



(9) 



It is interesting to derive the exponents A and n for 



Finally, in the limit 7 — >■ and a — ► 2, both Co and c\ go 
to zero, the distribution function is non-singular, n = 0, 
and the density profile after the adiabatic growth of the 
black hole has a 7 = 3/2 singularity, which is the result 
obtained in Ref. Q using an isothermal (non-singular) 
distribution function. 

A few comments on the results above. The derivation 
of the scaling in Eq. (||) verifies that the dark-matter 
spike is mainly formed by particles which, before the 
black-hole growth, were very cold and close to the Galac- 
tic center; i.e., particles with Ei — > <f>i(0). If the distribu- 
tion function is singular in that limit (i.e., there are many 
such particles in the initial state), the spike is hugely 
enhanced. As we have seen, this is possible even for a 
dark-matter-density profile that is non-singular towards 
the Galactic center (but in this case the spike is always 
milder than for 7 / 0). Note also that the enhancement 
in the spike due to the singularity in the distribution 
function is possible just if the black hole grows exactly 
at the center of the initial dark-matter system. If the 
black hole is formed away from the Galactic center, the 
induced dark-matter spike will be mild. For a system 
with an isotropic velocity-dispersion tensor, it will be of 
the same form we obtained above for a non-singular dis- 
tribution function. 
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In the previous paragraph, we found for the NFW pro- 
file a spike with A = 7/3. For this result to hold, it is 
crucial to postulate that the initial density profile has a 
1/r singularity all the way down to the Galactic center. 
This behavior is just an extrapolation from an N-body 
simulation with resolution not better than a kpc. As the 
simulations show that larger structures form from the 
merging of smaller substructures, such an extrapolation 
to scales smaller than 10 pc or so may not be trivial. 
More importantly, even if the 1/r scaling is indeed pre- 
served down to the inner few pc in N-body simulations, 
it is unlikely to survive the addition of baryons. If the 
density profile is assumed to be truncated in the inner- 
most region, we expect a shallower spike and a model- 
dependent suppression. For example, if the truncation of 
a singular profile is performed introducing a core radius, 
by replacing pi (r) oc 1/r 7 with pi(r) oc 1/[1 + (r/a) Q ] 7 / Q , 
the spike exponent A will depend just on a (its expression 
can be recovered from the analysis above substituting 
7 = and (3 = 7; note, however, that the normalization 
of the spike is set by the other parameters as well). 

Finally, the results in this Section show that the spike 
obtained by assuming a single-power-law density profile 
(~ 1 /r 7 for any r) with 7 — > will probably not provide 
a good approximation to the true spike. Most likely, the 
results in Ref. |l2|, where values of 7 as low as 10~ 6 are 
considered, will be modified assuming the more general 
double-power-law density profile of the form in Eq. (^J) . 

IV. INSTANTANEOUS GROWTH OF THE 
BLACK HOLE 

We now consider what happens when the black hole is 
inserted instantaneously to the center of the dark-matter 
distribution. This will provide some indication of what 
happens to the dark-matter distribution if the black- hole 
growth is too rapid to be approximated as adiabatic. We 
do this calculation under the assumption that all particles 
are on circular orbits. 

We first provide a qualitative argument to show that 
the sudden black-hole growth leads to a central dark- 
matter-density enhancement that is not as dramatic as 
that for adiabatic growth. If the black hole grows adia- 
batically, the initial large circular orbit becomes a final 
circular orbit of much smaller radius. Thus, the dark 
matter that was initially at a large radius spends all of 
its time at a much smaller radius. Now suppose the black 
hole appears suddenly and dominates the potential at 
these radii. Again, the angular momentum of the par- 
ticle is conserved, but its energy changes in such a way 
that it follows an elliptical Keplerian orbit, as shown in 
Fig. [2J in which its largest radius is the initial radius. 
Although the orbit reaches to smaller radii than the cir- 
cular orbit that arises from adiabatic growth, the velocity 
of the particle at these small radii is large. As a result, 
the particle on this orbit spends most of its time at radii 



initial orbit 




FIG. 2. Initial and final orbits for sudden and adiabatic 
black-hole growth. If the black-hole grows adiabatically, the 
initial large circular orbit becomes a final circular orbit of 
much smaller radius. If the black hole appears suddenly, the 
initial circular orbit becomes an elliptical orbit shown (assum- 
ing that the final potential at these radii is dominated by the 
black hole). 

larger than the radius of the final orbit from adiabatic 
growth. 

We now quantify these arguments. When the black 
hole is added, the velocity of each particle remains un- 
changed, but the gravitational potential is changed sud- 
denly. The final profile is obtained by following each 
particle along its subsequent orbit. The probability for 
a particle that at the appearance of the black hole was 
at galactocentric distance tq with a velocity vq to be be- 
tween radii r and r + dr is 



T(r ,v ) T(r Q ,v ) w,-(r, r , t/ ) ' 

where T is half of the period of the particle, and v r is the 
radial component of its velocity. 

If all particles are on circular orbits, then the calcu- 
lation is simplified, as the initial radial velocities are all 
zero. We thus need to integrate the probability that a 
particle that is observed at radius r in the final configu- 
ration came from an orbit at radius tq in the initial con- 
figuration. Doing so, we find that the final radial profile 
is given by: 

1 /" w 1 1 

Pf( r ) = — / dr orlpi{r )—--— r, (11) 

r J rmin T(r Q ) v r (r,r ) 

where 

v r = ^2 [Q f (r ) - $/(r)] + v*(r ) (1 - rg/H), (12) 



G 



with v c the circular velocity in the initial state. The 
extremes of integration r mul and r max are the mini- 
mum and maximum initial radii for particles contribut- 
ing to the final profile at the radius r. For such radii, 
v r {r, r^nin) = v r (r, r max ) — 0. From the expression for v r , 
it follows that r m j n = r; i.e., a particle on a circular orbit 
shrinks in general to a rosette orbit with apocenter equal 
to the radius of the initial orbit. 

We compute p f for the initial NFW profile intro- 
duced in Section III A . Again, we assume = 
$/(r) — GA/bhA; the result is shown as a dash-dot curve 
in Fig. |]. Although the black hole steepens the profile, 
the steepening is not nearly as dramatic as if the black 
hole is introduced adiabatically. 

The final profile has a r -4 / 3 singularity, a result that 
can be understood analytically. In the small-r limit, the 
potential is dominated by the black hole, and T(ro) <x 



3/2 



while v r becomes 



2GMbh 
rr 



Vo-r) 1 



2A/, 



BH 



V r 



For the NFW profile, the halo mass Mf 
within the radius rg is proportional to rg, 
computed exactly, and the final profile is 



(13) 



Vo) enclosed 
fmax can be 



Pf (r) ~ - 



dm 



ro) 

4n, (14) 



Vr§-(rr )/(fcri ax )-r/(fcr max ) r 4 /3 

as found numerically and shown in Fig. |j} 

Heuristically, our circular-orbit model should general- 
ize to any orbit with a finite peri-to-apocenter ratio. The 
apocenter distribution is unchanged after the implanting 
of the black hole, and the pericenter is determined by 
the angular-momentum distribution, which is a power 
law differing from that of the circular orbits by a scaling. 
So we expect the pericenter to move in according to the 
same power law as for circular orbits and thus produce a 
density power law as before except for normalization. 



V. MODELS WITH OFF-CENTERED BLACK 
HOLE 

The black hole, or at least its seed, might not form 
exactly at the center of the dark-matter halo due to de- 
viations from spherical symmetry in various galaxy for- 
mation processes. If the black hole forms off-center, then 
it will slowly sink via dynamical friction to the center of 
the Galactic potential well. We sketch in this Section the 
number of competing effects entering this problem. 
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FIG. 3. Predicted time for a seed black hole of mass 
Mf H = 10 4 M and 10 6 M to spiral in to the Galactic cen- 
ter from a radius r in an initially pure NFW or pure Moore 
halo profile (no baryons). It takes longer than a Hubble time 
for a 10 4 Mq black hole to spiral from r S> 100 pc, and for 
a 10 6 Mq black hole to spiral from r > 1 kpc. The spiral-in 
time is estimated with Chandrasekhar's formula for dynami- 
cal friction, depending slightly on the Coulomb factor and the 
eccentricity of the orbit of the black hole. Note the curves join 
together after a break at very small radius when the enclosed 
ambient mass becomes comparable to the mass of the infalling 
hole. Inside this radius, the hole would fall in quickly on a 
dynamical time scale, which is independent of M| H . 



A. Enough time to spiral in? 

The process of dynamical friction is not efficient 
enough to center the black hole regardless of its initial 
mass and of its initial distance from the center of mass. 
The spiral-in timescale for a very small black hole or for 
a black hole born very far from the center of the Galaxy, 
may be longer than the age of the Galaxy itself. For a 
black-hole seed of mass 10 4 M , the timescale to spiral 
in the potential well of the NFW profile we considered is 
roughly comparable to the Hubble time if its initial lo- 
cation is 100 pc away from the center of the dark-matter 
system. This is illustrated in Fig 0. Also shown is the 
case of a black- hole seed of mass 10 6 M Q , for which the 
timescale is of the order of 10 Gyr if it spirals in from 
1 kpc. The conclusion is that if, as plausible, the black- 
hole seed is generated through some violent process off- 
center, it will not get to the center until almost fully 
formed. 
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B. Merger and growth at the center 

What happens when the seed BH arrives at the cen- 
ter? There are two competing effects with the dark- 
matter distribution near the Galactic center: (a) it may 
be smoothed out by the time-dependent perturbation in- 
duced by the incoming black-hole; (b) it may be squeezed 
in by subsequent at-the-center growth of the black hole. 

Effect (a) has been investigated with N-body simula- 
tions in Ref. fl7|| . It was shown there that the back- 
reaction to the spiral-in of the black hole leads to the 
formation of a weak density cusp, with an r~°' 5 singu- 
larity, regardless of the initial density profile. More pre- 
cisely, the same features are obtained in simulations with 
the King model (which has a core radius) and with the 
Hernquist profile (which is singular at the Galactic center 
and given by Eq. (|) with (a, £,7) = (1,4,1)). Further- 
more, the size of the weak-cusp region is found to be 
correlated to the black-hole mass, with roughly the total 
mass of the region affected by the black hole equal to 
the black-hole mass. The proposed explanation for this 
result is as follows: The sinking of the black hole induces 
a heating process that leaves no particle in the system 
with energy below some minimum energy E . We saw 
above in the case of adiabatic growth at the center of the 
halo that the spike in the final density profile was due 
to initially cold particles. If such particles are removed 
from the system, neglecting the dependence on angular 
momentum in the distribution function, one finds that 
the final density profile goes at small radii like [[l7j 

Pf(r) = 4tt f f f (E f )y/2 [E f -*f(r)] 



47T 



ff(E f )J~2 [E f -* f (r)] 



(15) 



We implement this result in our picture, with a simple 
scheme. We suppose we have an initial dark-halo profile 
of the form in Eq. (|4|). Then we postulate the off-center 
formation of a black- hole seed of mass Mg H , which we 
keep as a free parameter. Since we are supposing that the 
black-hole seed appears off-center, regardless of whether 
it forms quickly or slowly, it induces a very mild dark- 
matter spike around it; we neglect this effect. We model 
then the modification to the dark-halo profile due to the 
sinking of the black-hole seed in the center of the Galaxy 
by replacing the initial density profile with the following: 



Ant(r) 



Ro 

Po 1 — 

r 



0.5 r 



1 + (Ro/r a ) 



1 + (r/r s ) 
1 + (i? /a) Q l (/3 ~ 7)/Q 



7-0.5 



1 + (r/a)» 



(16) 



which is essentially a modified Zhao model with a break 
at both r s and a. The new length scale r s we introduce 
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FIG. 4. Modification of an NFW dark-matter-density pro- 
file due to the off-center formation of a black-hole seed of mass 
Mjfjj, its spiral in the center of the dark-matter system and 
its adiabatic growth to the present-day mass of the black hole 
at the Galactic center. The cases for a few different values 
for the black-hole-seed mass are plotted. p C orc is the max- 
imum WIMP density above which WIMPs are depleted by 
pair annihilations. 

is determined by requiring that the mass of the dark halo 
within r s to be equal to Afg H . 

Moving onto Effect (b) , we model the accretion of mat- 
ter onto the black hole so that Mg H grows to the mass we 
now observe Mbh as a slow adiabatic process, and com- 
pute the final dark-matter density, treating the adiabatic 
contraction in the circular-orbit approximation described 
above. 

In Fig. H we pl ot pf for the initial NFW profile intro- 



duced in Section HI A and for a few choices for the mass 
of the black-hole seed. As can be seen, if the black hole 
that spirals in has a mass comparable to (more than a 
tenth of) the mass we observe now at the center of the 
Galaxy, the dark-matter density close to the black hole 
is, in this picture, lower than the density in the initial 
profile. This indicates that the number of dark-matter 
particles expelled from the inner Galaxy by the black- 
hole seed is larger than the number of those attracted 
later on by the adiabatic deepening of the potential well 
due to the increase in mass of the black hole. 

The adiabatic contraction process dominates for 
M| H < Mbh/10. For M| H < M BH /200, the end result 
is a dramatic enhancement in the dark-matter density 
around the black hole. Although the final density pro- 
file obtained with this procedure differs from what we 
found in case of adiabatic growth of the black hole at 
the center of the dark-matter system, the difference oc- 
curs only at densities higher than the maximum WIMP 
density /c C orc- Thus, if a black-hole seed with less than 
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FIG. 5. The same as in Fig. ^, but for a Moore et al. 
dark-matter-density profile. Also shown is the modification 
of the profile for sudden growth of the black hole at the cen- 
ter of the dark- matter system (dashed curve). 

I /200th the mass of the final black hole is let to grow adi- 
abatically at the center, then a spike that is essentially 
indistinguishable accordance with the GS result could be 
formed. The trouble is that such small seed cannot spiral 
in the center quickly enough, as shown by Fig. ||. So the 
conclusion is that if the GS spike is to grow, the black 
hole must form within the inner 50 pc of the dark-matter 
distribution: If it forms outside this radius with a large 
mass, it will destroy the cusp as it spirals in, and if it 
forms well outside this radius with a small mass, it will 
not have enough time to spiral in. 

Although we have considered so far just the NFW 
profile, the procedures we outlined can be applied to 
any density profile. As another example, we consider 
the dark-matter-density profile obtained in the high- 
resolution N-body simulation of Moore et al. jll). In 
this case, the initial dark-matter profile is more cuspy 
towards the Galactic center and is given by Eq. ([I]) with 
(a,/?, 7) = (1.5, 3, 1.5) (hereafter the Moore et al. pro- 
file) . Again we make a choice for the values of our galac- 
tocentric distance, the local halo density, and the length 
scale in agreement with available dynamical constraints: 
Ro = 8 kpc, po = 0.3 GeV cm -3 , and a = 28 kpc. 

Results in this case are shown in Fig. |^. The initial pro- 
file is shown as a dotted curve; the dashed curve shows 
the slight increase in density in case of instantaneous 
growth of the black hole according to the treatment in 
Section IV. The solid curves give the final profile in the 
case described above in this Section. We see again that 
for large values of the mass of the black-hole seed, the 
net effect is a decrease in the dark-matter density. For 
values Mg H < A/bh/50, the enhancement in the final 
density reproduces, for pf < p COIC , the result one gets 



just assuming the adiabatic growth at the center of the 
dark matter system. The timescales to re-center black 
hole seeds of 10 4 M Q and 10 6 M Q are again shown in 
Fig |^, and the same discussion we outlined for the NFW 
case applies here. 



VI. TAKING STARS INTO ACCOUNT 

So far we have used a schematic picture of the Galaxy 
which included just dark matter and the black hole, ne- 
glecting other baryonic components. More realistically, 
the stellar components may outweigh the dark matter in 
the inner part of the galaxy in the initial configuration 
before the black hole grows. We treat this problem here. 

In the standard sequence of galaxy formation, a 
roughly spherical overdensity, in which the baryons and 
dark matter are distributed homogeneously, undergoes 
gravitational collapse. The coUisionless matter then re- 
laxes to an NFW or Moore profile, as indicated by 
N-body simulations, while the baryonic matter shocks, 
heats, and then cools to a disk and/or bulge in the grav- 
itational potential well generated by both the dark and 
the baryonic matter. The radiative cooling and collapse 
of baryons is likely to be a slow process and hence it 
should be a good approximation to treat it as an adia- 
batic contraction of the baryons. 

We sketch this process in one example. We reconsider 
again the case of off-center formation discussed in the pre- 
vious Section, refining the definition of initial and final 
states. In the initial state we add the extra component 
due to baryons (same density profile as for dark matter 
but scaled down of a factor f^/SlA/)- In the final state we 
add the stellar bulge and disk. We pick standard forms 
for them, compatible with available observations in the 
outer part of the Galaxy and with the rotation curve 
close to the black hole (we do not reproduce the func- 
tional forms here as they do not enter critically in the 
discussion below). The spiral-in of the black- hole seed 
is treated in the same way as in Section [v| The adi- 
abatic contraction is performed using the circular-orbit 
approximation. We also replace the disk with a spherical 
mass distribution that produces the same rotation curve. 
This simplification allows us to still treat the problem as 
spherically symmetric and does not enter crucially in the 
results we focus on. 

Results are shown in Fig. |(| The adiabatic contrac- 
tion of the baryons in the bulge and the disk induces 
a significant enhancement of the dark-matter density at 
intermediate galactocentric distances (few pc up to few 
kpc). In the innermost region, however, the main role is 
still played by the black hole. 
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FIG. 6. The same as in Fig. ^, but taking into account the 
adiabatic contraction of baryons that form the bulge and the 
disc 

VII. CONCLUSIONS 

We have tried to determine the initial conditions re- 
quired for growth of this spike by considering several al- 
ternative formation scenarios. We have found that for 
the spike to arise, the majority of the black-hole mass 
must grow on a timescale long compared with 10 7 years 
and within 50 pc of the center of a cuspy dark-matter 
distribution. Moreover, the dark-matter particles that 
make up this spike must be initially in very cold orbits 
near the Galactic center. A spike essentially equivalent to 
the GS spike could in principle arise if a low-mass black 
hole spiraled in to the Galactic center and then grew in 
mass, but the spiral-in timescale for such low-mass black 
holes is larger than the Hubble time. 

If we understood precisely how the black hole in the 
Galactic center formed, we could say more definitively 
whether the GS spike arises or not. Unfortunately, how- 
ever, the details of black-hole formation are far from 
clear. The present-day growth rate for the Milky Way's 
black hole as well as the peak of AGN activity at redshifts 
z ~ 2 suggest that a good fraction of the black-hole mass 
was in place at early times. The most promising mech- 
anisms for black-hole production at high redshift do not 
seem to contain the ingredients required for the GS spike 
(see, e.g., Refs. JTsj). N-body simulations show that at 
that time, dark halos are forming via rapid and often vio- 
lent mergers. Moreover, hydrodynamic simulations show 
that when two galaxies collide, they can leave behind a 
cold dense baryonic gas that can then collapse to form a 
supermassive black hole. The rate for black-hole growth 
is estimated to be orders of magnitude larger than the 
current rate (see, e.g., Ref. ||l!|), although a timescale 
long compared with 10 7 years, the period for the orbits 



that eventually wind up in the spike, is still conceivable. 
However, even if the timescale is reasonable, the dark 
halos for each of the two initial galaxies pass through 
each other almost undisturbed and then later relax to 
the dark-matter profile for a single halo. The timescale 
for this relaxation should be even longer. Moreover, it 
is hard to believe that the cold orbits required for the 
GS spike would be found right at the center of the dark- 
matter distribution so soon after this very violent col- 
lision. The mass density in these inner regions will be 
dominated by the baryons, and baryons can shock, form 
molecular clouds, stars, star clusters, and supernova- 
driven winds. The dynamical effects of all these processes 
are heuristically expected to disrupt cold orbits. Finally, 
it is not clear why the centers of the dark-matter and gas 
distributions should be so coincident in this violent for- 
mation. (Although the black hole in our Galaxy coincides 
fairly well at the current epoch with what appears to be 
the dynamic center of the Galaxy, this is not necessarily 
the case in other galaxies, such as M31.) 

An alternative (although not mutually exclusive) pos- 
sibility is that a 10 3 — 10 5 M Q black hole could form at a 
redshift z ~ 30 when molecular-hydrogen cooling permits 
the collapse of the first baryonic objects. The remaining 
10 6 M of mass could then be accreted during the sub- 
sequent generations of mergers that eventually give rise 
to a galactic halo. Although it would not be unusual 
for the initial high-density peak in the primordial mass 
distribution that gave rise to this black-hole seed to be 
near the center of the resulting galactic halo, it would be 
somewhat surprising if it was so close to the center. 

In conclusion, we have determined the initial condi- 
tions for the growth of the GS spike and given several 
reasons to doubt they were manifest in our own Galaxy. 
Given these doubts, null searches for annihilation radia- 
tion from the Galactic center should not be interpreted 
as evidence against WIMP dark matter, nor against a 
cuspy halo. 

ACKNOWLEDGMENTS 

We thank A. Loeb for useful discussions. This 
work was supported in part at Caltech by NSF AST- 
0096023, NASA NAG5-8506, and DoE DE-FG03-92- 
ER40701, and at Trieste by the RTN project under grant 
HPRN-CT-2000-00152. H.S.Z. wishes to thank Institut 
d'Astrophysique de Paris and Sterrewacht Leiden for fi- 
nancial supports and Caltech for hospitality during his 
visit. 



[1] G. Jungman, M. Kamionkowski, and K. Griest, Phys. 
Rep. 267, 195 (1996). 



10 



L. Bergstrom, Rept. Prog. Phys. 63, 793 (2000 ) 



J. Kormendy and L. C. Ho, ^stro-ph/0003268 

A. Eckhart and R. Genzel, Nature 383, 415 (1996); 

B. M. Klein et al, Astrophys. J. 509, 678 (1998); R. 
Genzel et al., astro-ph/0001428. 
P. J. E. Peebles, Gen. Rel. Grav. 3, 63 (1972). 
P. Young, Astrophys. J. 242, 1232 (1980). 
J. R. Ipser and P. Sikivie, Phys. Rev. D 35, 3695 (1987). 
G. D. Quinlan, L. Hernquist, and S. Sigurdsson, Astro- 
phys. J. 440, 554 (1995). 

P. Gondolo and J. Silk, Phys. Rev. Lett. 83, 1719 (1999). 
J. F. Navarro, C. S. Frenk, and S. D. M. White, Astro- 
phys. J. 462, 563 (1996) [astro-ph/9508025]. 
B. Moore et al., Astrophys. J. 499, L5 (1998). 
P. Gondolo, hep-ph/0002226. 

G. Bertone, G. Sigl, and J. Silk, astro-ph/0011553; astro- 
ph/0101134. 

J. Magorrian and S. Tremaine, Mon. Not. R. Astron. Soc. 
309, 447 (1999). 

J. Binney and S. Tremaine, Galactic Dynamics (Prince- 
ton University Press, Princeton, 1987). 

H. S. Zhao, Mon. Not. R. Astron. Soc. 278, 488 (1996). 
T. Nakano and J. Makino, Astrophys. J. 525, L77 (1999). 
R. D. Blandford, astro-ph/9906025; M. J. Rees, astro- 
ph/9912346. 

[19] H. S. Zhao et al., in preparation. 



11 



